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I. INTRODUCTION 

Quantum entanglement is possibly the most intriguing 
property of states of composite quantum systems. It mani- 
fests itself in correlations of measurement outcomes that are 
stronger than attainable in any classical system. The renewed 
interest in a general theory of entanglement in recent years is 
largely due to the fact that entanglement is conceived as the 
key resource in protocols for quantum information process- 
ing. Initial investigations focused on the properties of bipartite 
entanglement of finite dimensional systems such as two-level 
systems. In fact, significant progress has been made, and our 
understanding of the entanglement of such systems is quite 
well developed [1]. A natural next step is the extension of 
these investigations to multi-partite systems. Unfortunately, 
the study of multi-partite entanglement suffers from a prolif- 
eration of different types of entanglement already in the pure 
state case [2], and even less is known about the mixed state 
case. For example, necessary and sufficient criteria for sepa- 
rability are still lacking. For other properties, such as distilla- 
bility, no efficient decision methods are known, and it is even 
difficult to find meaningful entanglement measures [3]. A di- 
rection that promises to lead to simpler structures is that of 
infinite dimensional subsystems, such as harmonic oscillators 
or light modes, which are commonly denoted as continuous- 
variable systems [4, 5]. Indeed, for continuous-variable sys- 
tems the situation becomes much more transparent if one re- 
stricts attention to Gaussian states (e.g. coherent, squeezed 
or thermal states) which are, in any case, the states that are 
readily experimentally accessible. 

Quite recently, it has been realized that it might be a very 
fruitful enterprise to apply the methods from the theory of en- 
tanglement not only to problems of quantum information sci- 
ence, but also to the study of quantum systems that are typi- 
cally regarded as belonging to statistical physics, systems that 
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consist of a large or infinite number of coupled subsystems 
[7, 8, 9, 10]. Examples of such systems are interacting spin 
systems, which, like most interacting systems, exhibit the nat- 
ural occurrence of entanglement, i.e., the ground state is gen- 
erally an entangled state [7, 8, 9]. It has, furthermore, been 
suspected that the study of the entanglement properties of such 
systems may shed light on the nature of the structure of clas- 
sical and quantum phase transitions [7, 9]. It has turned out, 
however, that the theoretical analysis of infinite spin chains is 
very complicated and only very rare examples can be solved 
analytically. Coupled harmonic oscillator systems allow for 
a much better mathematical description of their entanglement 
properties than spin systems. Physical realizations of such 
systems range from the vibrational degrees of freedom in lat- 
tices to the discrete version of free fields in quantum field the- 
ory. This motivates the approach that we have taken in this 
work, namely to investigate the entanglement structure of in- 
finitely extended harmonic oscillator systems. 

In this paper we study a special case, namely a set of har- 
monic oscillators arranged on a ring and furnished with a har- 
monic nearest-neighbor interaction, i.e., oscillators that are 
connected to each other via springs. The paper is organized 
as follows. In Section II we provide the basic mathematical 
tools that are employed in the analysis following in the re- 
maining sections. We then move on to derive a simple ana- 
lytical expression for the ground state energy of the harmonic 
oscillator systems. Our main interest is the computation of 
entanglement properties of the ground state of the chain. In 
Section III we derive a general formula for the logarithmic 
negativity [11, 12] which we employ as our measure of en- 
tanglement. In Section IV we present analytical results that 
concern the symmetrically bisected chain, that is the situa- 
tion where the chain is subdivided into two equal contiguous 
parts and the entanglement is calculated between those parts. 
We show how to construct a very simple lower bound on the 
log-negativity, in the form of a closed-form expression based 
on the coupling strengths; that is, no matrix calculations are 
necessary. Furthermore, for nearest-neighbor interaction, we 
show that the bound is sharp, i.e., gives the exact value of the 
log-negativity. Surprisingly, the value of the log-negativity in 
this case is independent of the chain length; in particular, it 
remains finite. We show in Section V that the problem is not 
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reducible to a four-oscillator picture, thereby demonstrating 
the non-triviality of the physical system. We then move on 
to Section VI, where we study general bisections of the chain 
numerically. We demonstrate that entanglement is maximized 
for the symmetrically bisected chain. Furthermore, and rather 
counterintuitively, for asymmetric bisections where one group 
of oscillators is very small, and especially when it consists of 
only one oscillator, we find that the entanglement decreases if 
the size of the other group is increased. We also demonstrate 
that for large numbers of oscillators the mean energy of the 
ground state and the value of the negativity are proportional 
and provide an interpretation for this result. In Section VII we 
discuss our results. We also provide an intuitive picture that 
allows to explain the results in the previous sections. 

Generally we have attempted to structure the sometimes 
somewhat involved mathematics in such a way, that the reader 
can skip it and extract the main physical results easily. We 
state at the beginning of each section what main result will be 
obtained and we state this result clearly, either in the form of 
a theorem or at the end of the section. 



II. COVARIANCE MATRIX FOR GAUSSIAN STATES OF 
THE HARMONIC CHAIN 

In this section we derive an expression for the covariance 
matrix of the ground state and of the thermal states of a set of 
harmonic oscillators that are coupled via a general interaction 
that is quadratic in the position operators (e.g. oscillators cou- 
pled by springs). As a byproduct we also give an expression 
for the energy of the ground state. 

Let us first consider the covariance matrix for the ground 
state of a single uncoupled harmonic oscillator. The Hamilto- 
nian is given by (we have adopted units where ft = 1) 

1 at mui 2 » 9 

H = — P 2 + X 2 . 

2m 2 

Denoting the quadrature operators as a column vector R, with 
Ri = X and R 2 = P, the Hamiltonian can be concisely 
rewritten as 

tt n>T ( mu 2 /2 \ 

H - R { l/(2m) ) K 

The covariance matrix 7 of a general state p is given by 

7M = Re Tr[p(R k - Tr[pR k ]){Ri - Tr\pRi])]. 

for 1 < k, I < 2. For p n the n-th eigenstate of the Hamilto- 
nian, p n = \n){n\, it is a straightforward exercise to calculate 
that 

We will only be interested in the ground state, p = |0)(0|, 
however, since this is the only eigenstate which is Gaussian. 

Passing to the harmonic chain consisting of n harmonic os- 
cillators, we will only consider interactions between the os- 
cillators due to a coupling between the different position op- 
erators. According to the (g,p)-convention we have adopted 



here, the vector R of quadrature operators is given by Rj — 

Xj and R n +j = Pj, for 1 < j < n. The Hamiltonian is then 
of the form 

fj r>T ( Vm^/2 \ 

H ~ R { K/(2m)) R ' 

where the n x n-matrix V contains the coupling coefficients. 
The Hamiltonian is thus written as a quadratic form in the 
quadrature operators; we will call the matrix corresponding 
to this form the Hamiltonian matrix (as opposed to H, the 
Hamiltonian operator). In the present case, the Hamiltonian 
matrix is a direct sum of the kinetic matrix l„/(2m) and the 
potential matrix Vmui 2 /2. 

In this paper, we will consider a harmonic chain "con- 
nected" end-to-end by a translationally invariant Hamiltonian. 
The y-matrix of the Hamiltonian is, therefore, a so-called 
circulant matrix [14]. This is a special case of a Toeplitz 
matrix because not only do we have Vj. k = Vj- k , but even 
Vj.k = v (j-k) mod n f° r 1 < j,k < n, due to the end-to-end 
connection. We can easily write the coefficients v k in terms 
of the coupling coefficients. For a nearest-neighbor coupling 
with "spring constant" K, the potential term of the Hamilto- 
nian reads 

n 2 

— 7T- Xl + K(X {k+1 ) mod „ - X k ) 2 . 

fc=l 1 

Therefore, we have 

v = 1 + «/(ma) 2 ), vi = -2K/{mLU 2 ). 

More generally, including fc-th nearest-neighbor couplings 
with spring constants K k , and defining 

2K k 
a-k = o ' 

we have 

v = 1 + 2(ai +a 2 + ...), 
Vj = —aj, for j > 0. 

The calculation of the corresponding covariance matrix can 
now proceed via a diagonalisation of the Hamiltonian ma- 
trix, which effectively results in a decoupling of oscillators. 
Since the commutation relations between the quadrature op- 
erators must be preserved, the diagonalisation must be based 
on a symplectic transformation S G Sp(2n, M). This means 
that we can only use equivalence transformations C ^ C — 
S T CS such that S T Y>S = S, where, in the (q,p) -convention, 
the symplectic matrix E is given by 

^ {-In 

This real skew-symmetric matrix incorporates the canoni- 
cal commutation relations between the canonical coordinates. 
Fortunately, because the kinetic matrix is a multiple of the 



3 



identity, the Hamiltonian matrix can be diagonalized by an 
orthogonal equivalence of the form 

C^C = {S® S) T C(S © S), 

where S is the real orthogonal n x n-matrix that diagonalizes 
the potential matrix V. It is readily checked that the resulting 
transformation is indeed a symplectic one. In fact, S® S is an 
element of the maximal compact subgroup of Sp{2n 1 R). 

So C is now diagonal and of the form C = {muj 2 /2)V 
l n / (2m), where V is the diagonal n x n-matrix with entries 
1 < j < n > the eigenvalues of V. The co variance matrix 
7' of the ground state of the transformed Hamiltonian consists 
therefore just of single-oscillator covariance matrices with pa- 
rameter u>j = uj^/rjj, 1 < j < n, and is diagonal itself, to wit, 

i = (7>7 P )/2, 

The covariance matrix 7 in the original coordinates is then 
obtained by transforming 7' back, 

7 = (s® sW(s®s) T 

= [(s 7 ^ T )e(s 7 ;s T )]/2 

= [(r 1/2 /H)ffi(^ 1/2 )]/2. 

To simplify the notation, we will henceforth set m = 1 and 
cj = 1. So we have a simple formula for the covariance matrix 
in terms of the potential matrix V, 

7 - (7z©7p)A 

7, = v- lj \ 

Using this same derivation, we can also easily find a for- 
mula for the energy of the ground state. We will need this re- 
sult in Section VI, where we will compare the log-negativity 
of a state to its energy. Indeed, the ground state energy 
of a single oscillator is Hlu/2. In the decoupled descrip- 
tion of the ground state of the chain, the oscillators have en- 
ergy fkd^/rjj/2, with r]j, 1 < j < n, being the eigenvalues 
of the potential matrix V. The total ground state energy is 
E = (fkj/2) X)"=i Denoting Hlu/2 by E Q , we therefore 
have 

E = E Q Tr[V 1/2 }. 

Finally, we turn to Gibbs states corresponding to some tem- 
perature T > 0, the states associated with the canonical en- 
semble, given by 

p(0) = exp(-/3tf)/Tr[exp(-/3tf)], 

where [3 = l/T. Again, one can obtain the covariance matrix 
7(/3) of the state p([3) in a convenient manner in the basis in 
which the Hamiltonian matrix is diagonal. The 2n x 2n di- 
agonal matrix j'((3) can be obtained using the virial theorem: 



the mean potential energy and the kinetic energy of a single 
oscillator are identical and half the mean energy of the system 
at inverse temperature [3. Using this procedure one obtains 

7'(/3) = (Yx(/?)®7;(/J))/2, 

= ;^( 1 + cxp(/4)-i )' 

In the convention where m = 1, u = 1, one gets 

7(/3) = (7,(/3)©7 P (/?))/2, 
7*08) = U- 1 /2(l„ + 2(exp(/3U 1 /2 ) _ ]ln) -i^ 

7p (/3) = U 1/2 (l„ + 2(exp(/3U 1 / 2 )-]l„)- 1 ), 

for the covariance matrix of a Gibbs state in the original 
canonical coordinates. 



III. GENERAL FORMULA FOR THE LOGARITHMIC 
NEGATIVITY 

In this section we derive a general formula for the loga- 
rithmic negativity of a Gaussian state of n coupled harmonic 
oscillators with respect to a bipartite split, given the covari- 
ance matrix 7 of the Gaussian state. This set may consist of 
all n oscillators or of a subset of m < n oscillators. The 
only restriction is that the covariance matrix must be a direct 
sum of a position part -f x and a momentum part j p , i.e., there 
must be no correlations between positions and momenta. The 
resulting formula can be found at the end of this section. 

Let mi and be the sizes of the two groups of oscilla- 
tors the entanglement between which we wish to calculate, 
and let m = mi + m 2 < n. From Section II, we know that 
the covariance matrix 7 of the ground state of the harmonic 
chain is given by 7 = (7^ © 7 P )/2, where j x = V~ x l 2 and 
7 p = V 1 / 2 . In order to calculate the entanglement between 
two disjoint groups of oscillators in this state, we need to con- 
sider the covariance matrix associated with the reduced state 
of the m oscillators of the two groups. This covariance ma- 
trix /i - from now on also referred to as reduced covariance 
matrix - is given by the 2m x 2m principal submatrix of 7 
that consists of those rows and colums of 7 that correspond 
to the canonical coordinates of either group 1 or group 2. If 
m = n, meaning that the whole set of n oscillators is consid- 
ered, this step is not necessary. The reduced covariance matrix 
/i is again of the form 

M = (Mz ©Mp)/ 2 . 

where both fi x and /i p are m x m-matrices. 

Taking the partial transpose of a covariance matrix corre- 
sponds to changing the sign of the momentum variables of 
the oscillators in the second group. This operation maps the 
covariance matrix [i to 

M r = PpP 
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with 



P P x © Pp , P x 1 m : 



P p is a m x m diagonal matrix. Specifically, the j-th diagonal 
element of P p is 1 or —1, depending on whether the oscillator 
on position 1 < j < m belongs to group 1 or 2, respectively. 

The logarithmic negativity [11, 12] of a state is defined as 
the logarithm of the trace norm of the partial transpose of the 
state. The negativity is an entanglement measure in the sense 
that it is a functional that is monotone under local quantum 
operators [12, 13]. To date it is the only feasible measure of 
entanglement for mixed Gaussian quantum states. The def- 
inition of the logarithmic negativity can be easily translated 
into an expression which does not involve the state itself, but 
rather the covariance matrix of the state: as the trace norm 
is unitarily invariant, one has the freedom to choose a basis 
for which the evaluation of the trace norm becomes partic- 
ularly simple. More specifically, one may make use of the 
Williamson normal form [15] for the partial transpose of the 
covariance matrix. The problem of evaluating the logarithmic 
negativity is then essentially reduced to a single-mode prob- 
lem. This procedure gives rise to the formula [12] 



N=- 



2m 

E 

k=l 



log 2 (min(l,2|A fe (z£-V r )l)), 



where Afc(iE V r )> 1 < k < 2m, are the eigenvalues of 
i£~V r . E is the symplectic matrix 







Since E -1 = — E, we have to calculate the spectrum of 
the matrix B = —iT,P/j,P, giving 2m real eigenvalues 
Afc(.B) of B. Then the logarithmic negativity equals N — 
— Y^k2i l°g2 min(l, 2|Afe(.B)|). This formula can be further 
simplified due to the direct sum structure of \i = {\x x © fj, p ) /2. 
Simplification of B yields 



B = - 
2 





fix 



-Ppflppp 




The eigenvalue equation of a block matrix of this form reads 

= A 



X 
Y 



which is equivalent to the coupled system of equations Xv = 
Am and Yu — Xv. Substituting one equation in the other 
yields XYu = X 2 u, hence the eigenvalues of the block 
matrix are plus and minus the square roots of the eigen- 
values of XY. In particular, the eigenvalues of B are 
±(Xj(li x PplJ,pPp/4:)) 1 ^ 2 , 1 < j < m. Because of the ±-sign, 
taking the absolute value of the eigenvalues has the effect of 
doubling the eigenvalue multiplicity. Hence, 

m 

N = -^log 2 min(l, \j(n x P P HpP p )), 

which is finally the resulting formula of the logarithmic nega- 
tivity in terms of the matrices \i x and ji p . 




n-2 



n/2+l 



FIG. 1 : The symmetrically bisected harmonic chain. The oscillators 
1 to n/2 form group 1, the oscillators n/2 + 1 to n form group 2. 



IV. THE SYMMETRICALLY BISECTED HARMONIC 
CHAIN 

In this section we present exact analytical results for the 
log-negativity in a chain of n harmonic oscillators with a 
translationally invariant coupling. Moreover, we shall be in- 
terested here in the most symmetric case of calculating the en- 
tanglement with respect to a symmetric bisection of the chain. 
That is, the number n of oscillators should be even and the 
oscillators in positions 1 to n/2 constitute group 1, the others 
group 2 (see Figure 1). Hence, in the notation of Section III 
H = 7- 

Using the result of Section III, we find that the logarith- 
mic negativity of a symmetrically bisected oscillator chain, of 
length n and with potential matrix V, is equal to 

n 

W = -]Tlog 2 (min(l,A,(Q))), 
i=i 



with 



Q = V~ 1/2 PV 1/2 P 

P = K/2®(~K/2)- 



A. Symmetry properties of Q 

We begin our analytical investigations by studying a more 
general object than Q, namely the matrix 

R = G-^PGP. 

Here, G is a general real circulant matrix that is symmetric 
under transposition. Since G is symmetric and circulant, it 
can be written in 2 x 2 block form as 



G 



G' G" 
G" G' 



5 



We will first show that R exhibits the same block structure. 
Define the n x n flip-matrix F = F n as 

Fi.j — $i,n-\-\ — j. 

To simplify the notation, we will mostly refrain from mention- 
ing the size n of F; the mathematical context should make it 
clear which n is being used. 

Lemma 1 The matrix R can be written in 2 x 2 block form as 
R _ i A B 



B A ■ 



Proof. Since G is circulant and symmetric, FGF = G, which 
is true for every symmetric Toeplitz matrix. Also, FPF = 
-P holds. Writing 



R 



A B 
C D 



symmetry demands that C — FBF and D = FAF. Further- 
more, both G and P are also invariant under the F n i 2 © F n i 2 
symmetry. Hence, R exhibits this symmetry too, i.e., (F © 
F)R(F ®F) = R. Thus, A and B are invariant under F. □ 



Matrices with this block structure can be brought in block 
diagonal form using a similarity transform: 

Lemma 2 Let S = (P + F)/y/2. Then 

s (b B a )s = {a + bf)®(a-bf). 

Proof. This follows by direct calculation and noting that 

S- 1 = S and 



In F 

F -U 



/V2. 



□ 



We now specialize the above results for R to the matrix 

Q = V- p PV p P, 

with V again a general real symmetric circulant matrix. The 
power p remains hitherto unspecified. Note that any power 
of a symmetric circulant matrix is again symmetric circulant. 
Again, the n x n-matrix V can be written in 2 x 2 block form 

V V"\ 
V" V ) ' 

where both V and V" F are Hermitian. By Lemma 1, Q can 
similarly be written as the block matrix 



Q' Q 1 
Q" Q 



The following Lemma is crucial for the rest of the calcula- 
tions. 



Lemma 3 With the previous notations, Q' + Q"F = (Q' — 

Q'^)- 1 and 

det{Q' + Q"F) = exp(-pTr[Flog(V)]), 
det(<9' - Q"F) = cxp(+pTr[Flog(F)]). 

For p in the interval < p < 1, the following also holds: if 
V"F > 0, then Q' + Q"F < l n/2 , and ifV'F < 0, then 
Qi - Q"F < l n/2 . 

Proof. Consider SQS. On one hand, we have, by Lemma 2, 

SQS = (Q' + Q"F) © {Q 1 - Q"F), 

and similarly, SVS = (V + V"F) © (V - V"F). Also, 
SPS = F, as a short calculation shows. On the other hand, 
we also have 

SQS = SV- p PV p PS 

= SV- p SSPSSV p SSPS 
= (SVS)- p SPS(SVS) p SPS, 

and therefore 

SQS = {{V + V"F)- p © (V - V"F)- p )F x 
((V + V" F) p © (V - V"F) P )F 

= {(V + V"F)- p © (V - V"F)- p ) x 
((V - V" F) p © (V + V"F) P ) 

= (V + V"F)- p (V - V"F) P © 
(V - V"F)- p (V + V"F) p . 

Identifying the blocks in the two expressions for SQS, we get 

Q' + Q"F = (V + V"F)- p (V - V"F) P , 
Q'-Q"F = (V -V"F)- p (V + V"F) P , 

so that Q' + Q"F is the inverse of Q' - Q"F and 

det(Q +Q F )-[ dct{v , + v „ F) ) ■ 

Furthermore, log V = S (log(V + V" F) ®\og(V -V" F))S, 
hence 

Tr[F\ogV] 

= Tr[SFS(log(V + V"F) © Iog(V" - V" F))] 
= Tr[P(Iog(V" + V"F) © Iog(V" - V"F))\ 
= Trpog(V" + V"F) - Iog(V" - V"F)\ 

dct{V + V"F) 
~ ° g det(V - V"F) ' 

This then yields 

dct(Q' + Q"F) = exp(-pTr[Flogy]), 
det{Q'-Q"F) = exp(+pTr[FlogV]). 

Considering the second assertion, if V" F > 0, then V' + 
V"F >V- V"F, and, by Lowner's theorem [14], 

(V + V"F) P > (V - V"F) P 
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for < p < 1. Hence, for any vector x ^ satisfying 
an equation (V - V"F)p X = X(V + V"F)p X , it follows 
that A must be less than or equal to 1 (to see this, take the 
inner product of both sides with the vector x). Rearranging 
the equation to 

(V + V"F)- p {V' - V'Ffx = Ax, 

which is just [Q 1 + Q" F)x = Ax, yields that the A for which 
such an x exists are precisely the eigenvalues of Q' + Q"F. 
Hence, under the condition V" F > 0, the eigenvalues of Q' + 
Q"F are less than or equal to 1. Similarly, if V"F < 0, we 
proceed in an identical way to show that the eigenvalues of 
Q' — Q"F are less than or equal to 1. □ 



B. A lower bound on the negativity 

We will now apply Lemma 3 to the case p = 1/2 and V 
being the potential matrix of the oscillator chain to obtain a 
lower bound on the logarithmic negativity: 

Theorem 1 The logarithmic negativity of the bisected oscil- 
lator chain of length n obeys 

N> |TY[Flog 2 (y)]|/2. 

If V"F is semidefinite (i.e., either positive or negative 
semidefinite), then equality holds. 

Proof. To calculate the negativity we need the eigenvalues of 
Q with p = 1/2 that are smaller than 1. By Lemma 2, the 
spectrum of Q is the union of the spectra of Q' + Q" F and 
of Q' — Q"F. By Lemma 3, for V matrices satisfying the 
V"F > condition, the eigenvalues of Q smaller than 1 are 
the eigenvalues of Q' + Q"F. Furthermore, 



Tr[\og(Q' + Q"F)} 



logdet(Q' + Q"F) 
-pTr[Flog(V)]. 



Setting p = 1/2 then gives N = Tr[Flog 2 (y)]/2. On the 
other hand, if V"F < 0, it is the eigenvalues of Q' - Q"F 
that we need to consider. Since Tr[log(Q' — Q"F)] = 
+pTr[Flog(V)], we find 

7V = -Tr[Flog 2 (V)]/2. 

For general V"F, we first note that the general formula for 
the negativity can be written as 

iV = -Tr[log 2 min(l„,Q)]. 

For commuting X and Y, mm(X, Y) is the elementwise min- 
imum in the eigenbasis of X (and Y). By Lemma 2, we then 
have N = -Tr[log 2 min(l„, (Q' + Q"F) (Q 1 - Q"F))] 
and this is also equal to N = — Tr[log 2 min(ll„/ 2 , Q' + 
Q"F)] - Tr[log 2 min(l n/2 ,0' - Q"F)}. From Lemma 3 
we also know that Q' + Q"F and Q' — Q"F are each other's 
inverse. Hence, 



and, because the two arguments of min commute, TV = 
-Tr[min(log 2 (Q' + Q"F),log 2 (Q' - Q"F))}. Finally, the 
trace of a minimum is smaller than or equal to the minimum 
of the traces, so that 

N > max(Tr[log 2 (Q' + Q"F)], Tr[log 2 (Q' - Q" F))]. 

Because the two arguments of max are each other's negative, 
the maximum amounts to taking the absolute value of, say, 
the first argument. Hence, N > |Tr[log 2 (Q' + Q"F)}\ = 
|Tr[i ;l log 2 (V r )]|/2, where the last equality follows from the 
first part of the proof. □ 



For the nearest-neighbor Hamiltonian, V is of the form 



V = 



/ «0 


Vi 










Vl \ 


Vl 


vo 


Vl 












Vl 


vo 
















Vl 












Vl 


vo 


Vl 


V Vl 










Vl 


vo / 



so 



V = 



( v vi 
vi v vi 

'•. '•. 



\ 







V"F 



V ••• 

(vi 




V o ... 



(1) 



vi vo vi 
vi v ) 

' \ 




vi I 



(2) 



N 



-Tr[log 2 min(Q' + Q"F, Q' - Q"F% 



As vi < 0, V"F is obviously a (negative) semidefinite ma- 
trix. Therefore, the nearest-neighbor Hamiltonian satisfies the 
equality condition of the Theorem, and the logarithmic neg- 
ativity of the bisected harmonic chain with nearest-neighbor 
Hamiltonian equals N = |Tr[Flog 2 (V)]|/2. 



C. An explicit formula in the coupling 

The bound of Theorem 1 is actually a very simple one be- 
cause we can give an explicit formula for |Tr[Flog 2 (V)]| in 
terms of the coupling coefficients ctj, as follows. 

Theorem 2 For a translationally invariant potential matrix V 
with coupling coefficients ol\, ct2, ■ ■ ■ ct m , 

|Tr[Flog 2 (V)]| = log 2 (l + 4(ai + a 3 + . . .)). 

Note, in this formula, the absence of the coefficients with even 
index. 
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Proof. The eigenvalue decomposition of a general circulant 
n x n matrix V is very simple to calculate. For convenience 
of notation, we use matrix indices starting from zero instead 
of 1. Let Vk.i = Vk-i, then V = WACl, with the kernel 
matrix of the discrete Fourier transform, 

Q k j = exp(fc^^)/Vn, 

with < k, I < n — 1. This matrix is unitary and symmetric. 
The eigenvalues A k are related to vi via a discrete Fourier 
transform according to 

A k = } cxp kl vi. 

For real symmetric V this gives 

2n 2tt 
Afc = v + 2vi cos(fc — ) + 2v 2 cos(2fc — ) + • • • . 

n n 

It is now a straightforward calculation to obtain an expression 
for Tr[F\og 2 (V)}. First, 

(QFftt) M 

= ^2 exp(jfc— )Sj !n -i-f exp(-j'l — )/n 
i,j'=o 



n-l 



J2eM\Jk-(n-l-j)l] — )/ 

' * 11 



n-l 



2exp([7(fc + 0-(«-l)0] — )/«■ 



All elements are zero except those for which k + lis an integer 
multiple of n, i.e., either fc = / = 0orfc + ^ = n: 



(QFQt) = V ex P ((0j + 0) — )/n = 1 

' u 



and 



n-l 



(OFOt)„_ M = Vexp((jn-(n-l)0— )/n 
ft 

27Ti 

= cxp(— i ). 

n 

In the calculation of Tr[F \og 2 (V)] we only need the non- 
zero diagonal elements of f2F£lt, which are the (0, 0) and the 
(n/2, n/2) elements. Hence 



Tr[Flog 2 00] 
(A ) + 1. 

v + 2v x + 2v 2 + ■ 



2ixi 

= log 2 ( A o) + log 2 (A n/2 ) exp((n/2) ) 

n 



= log 2 



Inserting the relations between the elements of V and the cou- 
pling coefficients ctj 

v = 1 + 2(ai +a 2 + ...), 
vj = —ctj, for j > 0, 



yields the stated formula. 



□ 



For the nearest-neighbor Hamiltonian, the only non-zero ctj 
coefficient is a = ct\, giving rise to the following simple ex- 
pression for the logarithmic negativity. 

Corollary 1 For the nearest-neighbor Hamiltonian with cou- 
pling coefficient a > 0, the logarithmic negativity of the bi- 
sected chain of length n is given by 

N= ilog 2 (l + 4a). 

It is remarkable indeed that the negativity is independent of n, 
the chain length. 



D. Other potential matrices 

To conclude this section, we will prove that any other circu- 
lant symmetric potential matrix does not satisfy the equality 
condition of Theorem 2 so that the negativity will in general 
be larger than the lower bound and, moreover, dependent on 
the size n of the chain. Consider first a Hamiltonian with a 
nearest-neighbor coupling of strength ai and a next-nearest- 
neighbor coupling of strength a 2 , where ai,a 2 > 0. The 
matrix V"F is then of the form 



V"F = 



( 



V 



-ax 
-a 2 



-ct 2 










-a 2 



-a 2 
-ai J 



The non-zero eigenvalues of this matrix are those of the sub- 
matrix 



-Ct\ 

-ct 2 



-a 2 




v — 2v\ + 2v 2 



As its determinant is negative, — a\, it is not a definite matrix, 
hence neither is V"F. 

Generally, a k-th neighbor coupling i.e., a coupling be- 
tween oscillators k places apart, yields a matrix V"F which 
is Hankel [14] and has — on two skew-diagonals. If there 
is a k such that a k and ctk+i are non-zero but ctk+ 2 = 0, then 
V"F contains a 2 x 2 principal submatrix of the form 

-a k -a k+ i 
-ctk+i 

which is again not definite. Hence, in that case, V"F is not 
semidefinite. Now, if one fixes the interactions and then let n 
grow (which is exactly the setting here), there will always be 
some point when V" F will exhibit a zero skew-diagonal and, 
hence, is not semidefinite. 
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o 
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FIG. 2: There exist no symplectic transformations that decouple all 
but four oscillators from each other in the case of the bisected har- 
monic chain with nearest-neighbor Hamiltonian. 
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V. INEQUIVALENCE TO A FOUR-OSCILLATOR 
PROBLEM 

At this point, one might be tempted to think that the inde- 
pendence of the log-negativity of the chain length n, in the 
case of nearest-neighbor interaction, is a consequence of the 
presumption that the bisected harmonic chain of length n > 4 
with nearest-neighbor interaction is in fact equivalent to a 
much simpler problem: there could be an appropriate choice 
of basis of the Hilbert spaces of system 1 and 2, corresponding 
to a symplectic transformation, such that, in effect, only those 
four oscillators that are adjacent to the split boundary would 
be in an entangled state. The other n/2 — 2 oscillators of each 
system would then be in pure product states, thereby not con- 
tributing to the logarithmic negativity. This would mean that 
one could locally disentangle all but four oscillators with local 
symplectic transformations (see Figure 2). 

If this indeed were the case, then symplectic transforma- 
tions Si,S2 € Sp(n, R) would exist such that 

7 = (S 1 ®S 2 ) T 1 '(S 1 ®S 2 ) 

7 ' = (i n _ 4 /2) e 712 e 712 e (i„_4/2) 

(note that we are using a quadrature ordering convention here 
that is different from the one used in the rest of this paper). 
Here, 712 and 712 are 4 x 4-covariance matrices associated 
with the oscillators 1 and n on the one hand and n/2 and 
n/2 + 1 on the other hand, and l„_4/2 is the covariance ma- 
trix of the pure product states of the remaining n/2 — 2 oscil- 
lators of system 1 and 2, respectively. If for any n such a basis 
change could be performed, leading to the same covariance 
matrices 712 and 712, then the invariance of the logarithmic 
negativity of the bisected chain - the statement of Corollary 1 
- would follow as a trivial consequence. We will briefly show, 
however, that this is not the case. 
Consider the eigenvalues of B' = 



B' = ffi' 



- 1 p((1„_ 4 /2) 712 © 712 © (1„- 4 /2))p. 



FIG. 3: Positive eigenvalues of Q — 1 versus a, for a chain of size 
n = 20. From Lemma 3 it follows that the eigenvalues of Q come in 
reciprocal pairs; hence, the plot shows that in this case all eigenvalues 
of Q are either larger than 1 or smaller than 1. For other chain sizes, 
the eigenvalues behave in a similar way. This shows that the sym- 
metrically bisected chain typically cannot be reduced to the system 
depicted in Figure 2. 



The spectrum of the corresponding matrix Q' that enters in 
the formula for the negativity can easily be evaluated using 
the procedure mentioned in Section III. It is given by 

U {Q') = {1, 1,91,92,93,94}, 

where qi , . . . , m > 0, and 1 appears n — 4 times. 

Now we can confront this result with the spectrum of the 
matrix Q of the harmonic chain as is. A simple numerical 
calculation yields the values depicted in Figure 3. Since the 
eigenvalues of Q come in reciprocal pairs, we only show the 
eigenvalues larger than 1; furthermore, we subtract 1 from 
them and show the result on a logarithmic scale (in order 
to clearly distinguish all eigenvalues). In the case depicted, 
n = 20, we see that 10 eigenvalues are larger than 1, for any 
value of the coupling constant a. Furthermore, the 10 remain- 
ing eigenvalues are all smaller than 1 . This means that, in fact, 
1 is not included in the spectrum of Q, which is completely at 
variance with the result for Q' of the purported reduced chain. 
Hence, we arrive at the statement that not even a single os- 
cillator can be exactly decoupled from all the others by the 
application on an appropriate local symplectic transformation. 

This analysis shows that the coupled bisected chain with 
nearest-neighbor interaction can not be reduced to a problem 
of only two pairs of interacting oscillators. In Section VII - 
equipped with further results from numerical investigations - 
we will discuss these findings and present an intuitive picture 
of the correlations present in the ground state of this system 
of coupled oscillators. 
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FIG. 4: Logarithmic negativity iV of a harmonic chain bisected in 
groups of size ni and 712. The interaction is nearest-neighbor with 
coupling a = 20. 



VI. GENERAL BISECTIONS 

In this section we turn towards more general problems, ex- 
hibiting less symmetry. As these problems are much more 
difficult to solve analytically, we basically have restricted our- 
selves to numerical calculations and we only give analyti- 
cal results for small subproblems, valid in some asymptotic 
regime only. 



A. Asymmetrical bisections 

In Figures 4 and 5 we show the results of a numerical 
calculation for asymmetrically bisected chains with nearest- 
neighbor coupling. That is, the groups of oscillators have 
sizes ri\ 7^ n 2 . From these figures a number of features 
are immediately obvious. The most striking feature is the 
"plateau" in the entanglement that is reached whenever both 
groups are sizeable enough (say ni,n 2 > 10, at least in the 
presented case for coupling strength a = 20). Of course, 
when ni = n 2 , being the "diagonal" of the plot, we recover 
the result of Section IV that the log-negativity is independent 
of n = ri\ + n 2 . From these figures we are led to conjec- 
ture that, in the case of nearest-neighbor coupling, the value 
of log-negativity for ri\ = n 2 is an upper bound on the values 
for ni 7^ n 2 (not to be confounded with the result of Theorem 
1, which says that this value is a lower bound for all symmet- 
ric bisections with general circulant couplings). Moreover, 
for general circulant couplings, we conjecture that an upper 
bound on 7V(ni, n 2 ) is given by linim^oo N(m 7 to). Another 
feature is that when, say, n\ is kept fixed the log-negativity 
decreases with n 2 from a given value of n 2 onwards. This 
phenomenon is seen most clearly with small n\, particularly 
for ni — 1. We will endeavour an intuitive explanation of 








FIG. 5: Same as Figure 4, but seen from a different viewpoint. 



N(n 1 ,n 2 )/N(oo, 00) 




FIG. 6: Effect of coupling strength a on the convergence of the 
log-negativity towards its maximal value. Group size is kept 
fixed at 20 and group size n\ is varied. Shown is the ratio 
N(ni,ri2)/N(oo, 00). The different curves are for various values 
of q. One clearly sees that for small couplings the limit value is 
reached much faster. 



these features below. In conclusion, we conjecture that, again 
for general circulant couplings, lim„ 2 _^ 00 N(l, n 2 ) is a lower 
bound on N(m, n 2 ). 

From Figure 6 we can see that the convergence of TV to- 
wards its plateau value N(oo, 00) depends on the strength 
of the coupling a. For higher values, convergence is slower. 
What cannot be seen from this figure is that the actual plateau 
value is larger as well. 
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B. Entanglement versus energy 

It is interesting to compare the entanglement present in 
the chain ground state with its energy. We consider nearest- 
neighbor interaction only. We have shown in Section II that 
the ground state energy equals (fkj/2)Tr[V 1 / 2 ]. For zero 
coupling (V = 1„) this gives just n times the single-oscillator 
ground state energy E = Hui/2, as expected. For large cou- 
plings a, we show that the ground state energy is of the order 
of %JoluEq. 

From the proof of Theorem 2, we have that the eigenvalues 
Afe of V are given by 

Afe = vo + 2«i cos(fc27r/n) + 2v 2 cos(2fc27r/n) + . . . . 

The energy in terms of these eigenvalues is 2^fc=o ^l/ 2 - F° r 
large values of n, we can replace the discrete sum over k by 
an integral in x — 2nk jn. For nearest-neighbor coupling, this 
yields: 

E w E (n/n) [ dx(vo -\- 2v\ cos 2x) ^/ 2 
Jo 

,tt/2 

= 2E (n/ir) / dx((v + 2«i) - 4«i sin 2 x) 1/2 
Jo 

,tt/2 

= 2E (n/ir) / + 4asin 2 x) 1 / 2 

Jo 

,tt/2 

= 2E {n/-K)s/a I dx(4sin 2 (a;) + l/a) 1/2 . 
Jo 

In the limit of a tending to infinity, the latter integral tends 
to 2 J^ 2 dx sin(x) = 2, so that indeed 

4 

7T 

Recalling the exact formula for the log-negativity in the sym- 
metrically bisected case, we have that the negativity (not the 
logarithmic one) is \/l + 4a. We thus find that the negativ- 
ity is approximately proportional to the mean energy per os- 
cillator. The exact values, calculated numerically, have been 
plotted in Figure 7. For a = 0, the curve obviously goes 
through the point with mean energy equal to Eq and negativ- 
ity equal to 1 . For a going to infinity, the mean energy goes 
to (2/n)E a = 0.63662^0 times the negativity. 

C. Non-contiguous groups 

From the above, one would get the impression that the mean 
energy gives a general upper bound on the amount of entan- 
glement in the system (apart from a numerical factor). This 
is certainly not the case, because, until now, we only have in- 
vestigated the cases where the two groups of oscillators were 
contiguous. In the following paragraph we look into the en- 
tanglement between non-contiguous groups. Specifically, we 
look at the extreme case of entanglement between the group 
of even oscillators and the group of odd ones. As can be seen 



Energy per oscillator per unit of negativity 




1 23456789 10 11 



N 

FIG. 7: Energy per oscillator (in units of hto/2) per unit of negativity 
(not logarithmic) in function of the negativity, for the case of con- 
tiguous groups of large enough size (so that the entanglement plateau 
in Figure 5 is reached). The interaction is nearest-neighbor and the 
coupling is implicitly present as a parameter. The dashed line depicts 
the limiting value for infinitely strong coupling. Here, the number n 
of oscillators is taken to be 20. However, the results become inde- 
pendent of n for n large enough: for n = 20, the limiting value is 
0.63531, while for infinite n the exact result is 2/7T = 0.63662. 

from Figure 8, numerical calculations already show that the 
log-negativity tends to a constant times n, the chain length. 
Therefore, in this case, the log-negativity can grow indefi- 
nitely large even when the mean energy is kept fixed. In view 
of this, it would be more correct to say that there are two con- 
tributions to the entanglement: one is the mean energy, which 
is directly related to the coupling strengths, and the second 
is the surface area of the boundary between the two groups of 
oscillators, which in the 1 -dimensional case is just the number 
of points where the two groups "touch" each other. We will 
return to this issue in Section VII. 

The validity of the purported linear relationship can be 
shown analytically in a rather simple way, yielding as a by- 
product an expression for the proportionality constant. First 
of all, the diagonal elements of the P matrix for this configu- 
ration are 1 for odd index values, and -1 for even index values. 
The Fourier transform of P, that is: QPCl^ (see proof of The- 
orem 2), is equal to 

( K/2 \ 

V K/2 J ' 

as is easily checked. We already have calculated the Fourier 
transform of the V matrix, which again can be inferred from 
the proof of Theorem 2. It is given by A = flVrt; here, 
A is a diagonal matrix with diagonal elements A^ = v + 
2v\ cos(2/c7r/n) + 2w 2 cos(4fc7r/n) + •••, 0<fc<n — 1. 
Inserting this in the expression for the Q-matrix gives: 
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N between even and odd group vs.w 



N 




FIG. 8: Entanglement between the group of even oscillators and the 
group of odd oscillators, in function of the chain length n (even n 
only). Interaction is again nearest-neighbor with coupling a = 20. 
The log-negativity is seen to quickly converge to a constant times n. 
The value of the constant depends on a and the relationship is shown 
in Figure 9. 



If we write A in 2 x 2 block form as 



A' 
A" 



the spectrum of Q is the union of the spectrum of 
A ,-i/2 A „i/2 and Qf A //-i/2 A /i/2_ Worked . out5 gives 

the eigenvalues (A k+n/2 / Ak) 1/2 and (A fc /A fe+ „ /2 ) 1/2 , for 
< k < n/2 — 1. Using the inherent symmetry that 
A„_fc = Afe, the eigenvalues of Q are 



A* 



Vi/2-fc 



±1/2 



The formula for the log-negativity obtained in Section III can 
be reformulated as minus the sum of the negative eigenvalues 
of log 2 Q. In the present case we get as log-negativity 



,./2 



N : 



\ J2 1 lo §2 



A fe 



A, 



k=0 -n/2-k 

For the nearest-neighbor Hamiltonian, this simplifies to 



N 



n/2 
fe=0 



n/4 

E 

k=0 



log. 



v + 2vi cos(27rfc/n) 



log 2 



2 v — 2v\ cos(27rfc/n) 

l + 2a(l + cos(27rfc/n)) 
1 + 2a(l - cos(27rfc/n)) 



for n that are multiples of 4. For large n, we can replace the 
discrete sum by an integral, 



2tt 



* /2 , , l + 2a(l 
ax log 2 



cos (a;)) 



1 + 2a(l - cos^))' 



FIG. 9: Relationship between the constant factor c appearing in the 
asymptotic formula for the log-negativity N — cn in the even-odd 
setting (as in Figure 8) and the coupling constant a. 



which indeed proves that, for large n, the log-negativity is a 
linear function of n. The integral itself cannot be brought in 
closed form. In Figure 9, we show the result of numerical 
calculations giving the asymptotic value of N/n versus a. 



D. Effect of group separation 

In the following paragraph, we give some results for con- 
tiguous groups that do not comprise the whole chain. In Fig- 
ure 10 we consider a fixed chain of n = 40 oscillators and 
look at the entanglement between two equally sized contigu- 
ous groups, in function of the group size and the separation 
between them. We define the separation as the number of os- 
cillators in the smallest gap between the groups; since we are 
dealing with a ring, there are two gaps between the groups. 
Note that the log-negativity is plotted on a logarithmic scale. 
There are two main features in this figure. The first and least 
unexpected feature is that the entanglement decreases more or 
less exponentially with the separation. We believe that this is 
quite natural in view of the fact that the coupling between the 
groups also decreases with the distance. 

The more remarkable feature is that for small groups, the 
entanglement quickly becomes zero altogether, as measured 
by the logarithmic negativity. Bound entanglement [6, 17] is 
of course not detected by this measure of entanglement, and it 
would be an interesting enterprise in its own right to study the 
structure of bound entanglement present in coupled oscillator 
systems. We will leave this, however, for future investigations. 
From now on, the term that no entanglement is present will be 
used synonymically with the statement that the logarithmic 
negativity vanishes. 

For groups of size 1, the log-negativity is zero already at 
separation 1 (the gap consists of one oscillator). For groups 
of size 2, there is still entanglement at separation 1, but none 
at separation 2. The larger the groups, the larger the maxi- 
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N vs. group separation and group size s 
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FIG. 10: Log-negativity for two contiguous groups that do not com- 
prise the whole chain. The chain consists of 40 oscillators, coupling 
is nearest-neighbor with coupling strength a = 20. Shown is the 
log-negativity, displayed on a logarithmic scale, versus the separa- 
tion between the groups, i.e., the number of oscillator positions be- 
tween them. The different curves are for different group sizes s (both 
groups are taken to be equal in size). The curve for group size 1 is 
not visible because it is a single point: the log-negativity between 2 
oscillators turns out to be whenever their separation is larger than 
0. 



mal separation for which there is still entanglement can be. 
One could try to interpret this by saying that there is a kind 
of threshold value below which entanglement drops to zero. 
However, this is more a reformulation of the results than an 
explanation, because it sheds no light on why this supposed 
threshold should depend on the group size. 

To really explain what is happening, we need to take a 
closer look at the exact calculations. Consider first two groups 
of oscillators of size 2 and with separation 1 ; that is, group 1 
is at positions 1 and 2, group 2 at 4 and 5. The Q-matrix of 
this configuration (with n — 40 and a = 20) has eigenvalues 
2.063, 1.1339, 1.0938 and 0.88361. As one of the eigenval- 
ues is smaller than 1, there is, indeed, entanglement present. 
We might be led to think that this entanglement is the cumula- 
tive result of the entanglement between the different oscillator 
pairs, (1,4), (1,5), (2,4) and (2,5), but this is not true, because 
these pairs are not entangled themselves: their separation is 
larger than 0. What is happening here is that the eigenvalues 
of the Q-matrix belonging to pair (1,4), say, are 1.8065 and 
1.1724, which are both larger than 1 and, therefore, do not 
count in the entanglement figure. 

The resolution of this strange behaviour in terms of the sep- 
aration is that the mere fact alone of having correlations be- 
tween the groups (eigenvalues of Q different from 1) is not 
enough to have entanglement. The correlations must be of 
special nature, namely: the eigenvalues of Q must be smaller 
than 1 . One could say that larger groups can more easily ex- 
hibit entanglement; their Q matrix has a larger dimension and, 




FIG. 1 1 : Classical correlations in a chain consisting of 40 oscillators; 
coupling is nearest-neighbor with coupling strength a = 20. Shown 
is the quantity (XiXj), displayed on a logarithmic scale, versus the 
second oscillator index j. 



hence, more eigenvalues, so that there are more opportunities 
for having at least one eigenvalue smaller than 1 . 

In this respect, it is interesting also to have a look at the 
classical correlations in the chain, i.e. the expectation values 
(XjX k ). From the treatment in Section II we immediately 
see that these correlations are given by the elements of the 
matrix j x /2 = V~ 1 / 2 /2. In the case of circulant symmetry, 
we only have to consider the first row of the matrix, giving 
the correlations between the first oscillator and any other one. 
Figure 11 shows these classical correlations for the system 
considered in Figure 10 (n = 40, a = 20). As could be 
expected, these correlations decrease exponentially with the 
oscillator distance and, furthermore, never vanish completely. 



E. Thermal state 

To conclude this section, we consider a thermal state in- 
stead of the ground state. The calculations are exactly the 
same in both cases, apart from the fact that in the covariance 
matrix there is an additional factor l„+2(exp(/3-\/V 7 )- 1„) _1 
to the 7 P and j x blocks (T = 1//3). The results are shown in 
Figure 12. One sees that for small temperatures the negativity 
is equal to the ground state negativity, and from some value 
onwards it starts to decrease more or less linearly with T until 
there is no (free) entanglement at all anymore. 



VII. DISCUSSION 

The numerical results we have obtained in Section VI can 
be interpreted in a qualitative way, by means of two rules-of- 
thumb. These rules are not to be interpreted as strict math- 
ematical statements; for that, we already have the exact for- 
mulas. The importance of the two rules is that they allow to 
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FIG. 12: Log-negativity of a thermal state with temperature T versus 
T and chain size n. Symmetrically bisected chain, nearest-neighbor 
interaction with coupling a = 20. 

reason about the dependence of entanglement on various fac- 
tors, like group size, coupling strengths and group geometry. 

The first rule is that, due to the coupling between the oscil- 
lators, the system exhibits inter-oscillator correlations which 
are decreasing with distance. This is a fairly natural state- 
ment, in view of the fact that the couplings between the oscil- 
lators are short-range as well. In a more mathematical way, 
one could consider the matrix j x = y~ 1//2 , whose elements 
are the classical correlations (XjXk). The j-th row describes 
the correlations between the j-th oscillator and all other os- 
cillators. The correlations can thus, in a figurative way, be 
subdivided into packets, one packet for every row in the cor- 
relation matrix. For chains with a circulant potential matrix V, 
it is self-evident that j x is also circulant so that the correlation 
packets all have an identical shape (see Figure 13). 

The second rule is that the entanglement between two 
groups of oscillators depends on the total amount of correla- 
tion between the groups. Again, this rule looks fairly innocu- 
ous and even trivial. However, combining the two rules read- 
ily shows why, in the case of contiguous groups, the entangle- 
ment in function of the group sizes should reach a plateau. In- 
deed, even while the total amount of correlation grows, more 
or less, linearly with the chain size, this has very little impact 
on the entanglement between the groups because it are only 
the correlation packets that straddle the group boundaries that 
enter in the bipartite entanglement figure. For large groups, 
most of the correlation packets describe correlations within 
the groups. What is important is the amount of correlations 
between the groups, and this quantity is virtually independent 
on the group size, provided the groups are so large they can 
accomodate most of the packets within their boundaries. 

We must stress, however, that these two rules are of a qual- 
itative nature. As noted already in Section VI, in the discus- 
sion of the dependence of entanglement on group separation, 
having correlations between the groups alone is not enough 




FIG. 13: Schematic drawing of the inter-oscillator correlation pack- 
ets, i.e. the rows of the correlation matrix *y x . The line thickness 
indicates the amount by which the correlation packet is involved in 
the entanglement between the groups, i.e., how much it is shared by 
the two groups. 



for having entanglement. The correlations must be such that 
the Q matrix has at least one eigenvalue smaller than 1 . The 
bottom line is in any case that one must go through the exact 
calculations to see whether or not there is entanglement. 

Another effect that can be accounted for is the dependence 
of the log-negativity on the group size if at least one group is 
very small. If both groups are very small, say 1 oscillator both, 
then the packets are so wide they wind up along the chain and, 
therefore, cross every group more than once, adding to the en- 
tanglement figure a number of times. If one of the groups is 
kept fixed, and the other is made larger, the winding number 
of the packets decreases and so does the amount by which the 
packet is shared by the groups. This effect could explain the 
decrease of entanglement with growing n 2 . At this point, the 
qualitative reasoning again breaks down, however, since the 
reduction of the amount of sharing per packet is counteracted 
by an increase in the number of packets. To show that the bal- 
ance is still in favour of an entanglement decrease, once again 
one really needs to go through the exact calculations; this is 
what we have done in Section IV. Nevertheless, the qualitative 
reasoning has the virtue that it shows what the main ingredi- 
ents are. Furthermore, it immediately leads to the conjecture 
that the effect of decreasing entanglement would not occur in 
a chain that is not connected end-to-end, since no winding oc- 
curs there. 

A numerical experiment immediately showed that this is 
exactly what happens, as witnessed by Figure 14. One has 
to be careful, though, about how one "opens" the chain. To 
clearly show the disappearance of the winding effect, one has 
to make sure that opening the chain does not introduce side- 
effects. Particularly, the oscillators at both ends should still 
"see" the same springs as before opening the chain. One can 
take care of this by connecting the ends of the chain to two 
additional oscillators that are kept in a zero-energy state (i.e. 
with zero X-variance; hence, they must be oscillators with in- 
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FIG. 14: Same as Figure 5, but with a Hamiltonian that does not con- 
nect the harmonic chain end-to-end. The two groups are, therefore, 
connected only at one point (in the middle of the chain), and this ex- 
plains why the log-negativity is only about one half of the value it 
had with end-to-end connection. Furthermore, for small m we now 
see an increase with ri2 instead of a decrease, which seems to imply 
that the counterintuitive behaviour on the ring is actually a winding 
effect (see text). 

finite mass). At the level of the potential matrix, this means 
that the diagonal elements V\ i and V n n are still 1 + 2a, al- 
though the elements Vi )Tl and V n _\ are being set to zero. Not- 
ing the analogy between harmonic chains and transmission 
lines, we call this special connection process the termination 
of the harmonic chain. In transmission line theory, correct ter- 



mination of a line (using appropriately matched impedances) 
is necessary to avoid signal reflections at the ends of the line. 
We believe that analogous reflection effects could be exhib- 
ited by non-terminated harmonic chains, but leave the investi- 
gation of this boundary phenomenon to future work. 

Finally, for non-contiguous groups, and, specifically, for 
the entanglement between even and odd oscillators, the two 
rules-of-thumb correctly predict that the entanglement keeps 
increasing with growing chain length. Indeed, if n grows, 
then the "boundary area" between the even and odd group 
also grows (linearly with n), in contrast with the contiguous 
groups, whose boundary area is fixed (1 for the open chain, 2 
for the closed chain). Hence, the amount of correlations strad- 
dling the boundary should grow too. The exact calculation 
confirms this effect and shows a linear relationship between 
log-negativity and chain size. It would be interesting to inves- 
tigate what happens in three-dimensional oscillator arrange- 
ments with couplings decreasing with distance. We believe 
that a similar relation will show up between entanglement be- 
tween two groups and the area of the boundary between the 
groups. We leave this issue, however, for future investiga- 
tions. 
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